THE MAGNUS EXPANSION, TREES AND 
KNUTH'S ROTATION CORRESPONDENCE 

KURUSCH EBRAHIMI-FARD AND DOMINIQUE MANCHON 

Abstract. W. Magnus introduced a particular differential equation characterizing the logarithm of the solution of 
*v/ linear initial value problems for linear operators. The recursive solution of this differential equation leads to a 

peculiar Lie series, which is known as Magnus expansion, and involves Bernoulli numbers, iterated Lie brackets 
and integrals. This paper aims at obtaining further insights into the fine structure of the Magnus expansion. 
£\J ■ By using basic combinatorics on planar rooted trees we prove a closed formula for the Magnus expansion in 

the context of free dendriform algebra. From this, by using a well-known dendriform algebra structure on the 
vector space generated by the disjoint union of the symmetric groups, we derive the Mielnik-Plebahski-Strichartz 
formula for the continuous Baker-Campbell-Hausdorff series. 
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1. Introduction 

During the last decade, some suiprising convergence of different areas of mathematical sciences has oc- 
curred. The seemingly separate fields: numerical integration methods, Lyons' rough path theory E3l l36l . 
Ecalle's mould calculus Ifl2l . and Connes' noncommutative geometry |[T6l . share a common algebraic formal- 
ism where algebraic structures on trees and its underlying combinatorics are central. 

Let us give two closely related examples coming from the theory of numerical integration. First we men- 
tion the pioneering work of J. Butcher on an algebraic theory of integration methods in the 1960s and 1970s 
CJEH. Butcher's S-series can be seen as a generalization of Taylor series, in which rooted trees naturally 
appear. They uniquely represent elementary differentials, as Cayley noticed already in his classical 1857 paper 
llin . Motivated by the problem of extending Butcher's work to the construction of generalized Runge-Kutta 
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methods for integration of differential equations evolving on Lie groups, Munthe-Kaas introduced in RU1 |4T1 
the notion of Lie-Butcher series. In applications of both Butcher and Lie-Butcher series, algebraic and related 
combinatorial structures on rooted trees (in non-planar and planar version, respectively,) play an essential role. 
Since then, identifying genuine algebraic structures became a useful part of the theory of numerical integration 
methods, see e.g. |[T4llT5ll34ll35ll43l . The other example evolved out of W. Magnus' seminal 1954 paper 
(371, where the author introduced a particular differential equation characterizing the exponential solution of 
differential equations for a linear operator in terms of a Lie series. The latter is known as Magnus expansion 
and has become a well-known tool in the solution and approximation theory of linear initial value problems 
(5l|42j|46l. Iserles and collaborators E71 1281 |29l were the first to use planar tree structures in an intriguing 
way to study the Magnus expansion in the context of numerical integration. 

Recently it became clear that in general most of the combinatorial structures on trees can be traced back 
to the fact that free pre-Lie and dendriform algebras are naturally described in terms of rooted trees and pla- 
nar binaiy trees, respectively irT3ll3Tll32l . Indeed, trees provide genuine examples for combinatorial objects 
spanning connected graded locally finite-dimensional vector spaces, on which rich and various extra alge- 
braic patterns are given in explicit terms by so-called tree grafting operations. These algebraic structures are 
summarized by the notion of combinatorial Hopf algebra, see e.g. IT2611 . 

As it turns out, both Butcher's B-series and Magnus' expansion are most naturally described in terms of 
pre-Lie algebras. Indeed, Chapoton [|9| was the first to study B-series as genuine formal series of rooted trees 
in the free pre-Lie algebra in one generator. In fact, later it became clear that also Magnus' series is contained, 
though in more disguised form, in his 2002 preprint. See also Murua's work P3l for a link of Magnus' series 
to the Butcher-Connes-Kreimer Hopf algebra of rooted trees, and the underlying pre-Lie algebra. Based on 
the work by Iserles et al., we studied in |[TFl[T9l the Magnus expansion in the light of its underlying pre-Lie 
structure, using the corresponding dendriform algebra. Later we realized that Agrachev and Gamkrelidze [1] 
wrote down the Magnus expansion in the general context of chronological algebras, which is another name 
for pre-Lie algebras, as early as 1981. The interpretation as a logarithm, however, necessitates the dendriform 
structure as well. 

This paper is a continuation of our work in |TT8~I . It aims at obtaining further insights into the fine structure of 
the Magnus expansion by using basic combinatorial methods steaming from the description of free dendriform 
algebra in terms of planar binary trees. This allows us to present a closed formula for the Magnus expansion. 
A dendriform algebra structure on the linear span of the symmetric groups 11331 allows us to derive the contin- 
uous Baker-Campbell-Hausdorff series. In other words, we recover the known Mielnik-Plebahski-Strichartz 
formula for the (classical) Magnus expansion. Using fairly elementary tools, our result may be seen as a pedes- 
trian approach to parts of Chapoton's et al. irTOllTTI and Thibon's et al. IT221 work on the Magnus series, where 
descent algebra, operads, and the theory of non-commutative symmetric functions play a dominant role. 

The paper is organized as follows: in Section |2] we introduce the required structures, both combinatorial 
(planar binary trees, planar rooted trees and Knuth's rotation correspondence) and algebraic (Rota-Baxter, 
dendriform and pre-Lie algebras). We show how the dendriform algebra structure on the vector space spanned 
by planar binary trees together with the associative product is transported through Knuth's rotation correspon- 
dence to the vector space spanned by planar rooted trees. We describe the classical Magnus expansion in 
Section [3] writing it in terms of the pre-Lie product given on the space of locally integrable matrix- valued 
numerical functions by: 



X 

Jo 



f(u) du, g(s) 



(f>g)(s) 
The corresponding dendriform products are given by: 

(f > g)(s) := U /(«) du\ g{s), (f < g)(s) := f{s) U g(u) du) . 

In Section 0] we look for a closed formula (Theorem |U in the completion of the free unital dendriform algebra, 
giving the logarithm of the solution X of the linear dendriform equation: 

X = 1 + a < X. 
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A ^-analog of this result, relying on the abstract notion of descents for planar binary trees can be found e.g. in 
a recent paper by Chapoton iflOl . In our approach however, we refrained from using tools from the theory 
of non-commutative symmetric functions B2l[T7l . and derived the formula directly by invoking elementary 
combinatorial methods. Through the rotation correspondence the planar rooted tree picture reveals to be par- 
ticularly useful here, since the number of descents of a planar binaiy tree simply corresponds to the number of 
leaves of its planar rooted image, excluding the leftmost leaf. The formula in the free setting yields of course a 
similar formula in any complete filtered dendriform algebra (Corollary [6]). In the last section we recover from 
Theorem \4\ and Corollary [6] the well-known Mielnik-Plebahski-Strichartz formula ||42] 1461 for the classical 
Magnus expansion, also know as continuous Baker-Campbell-Hausdorff series. For that purpose we use a 
dendriform algebra structure on the direct sum ® >j k[S n ], where S n is the permutation group of n letters 
ETl 1331 . and the fact that the notion of descent for a planar binary tree matches well with the well-known 
corresponding notion for a permutation. 

Acknowledgements: We thank H. Munthe-Kaas and A. Lundervold for discussions and remarks. The first 
author is supported by a Ramon y Cajal research grant from the Spanish government. Both authors were sup- 
ported by the CNRS (GDR Renormalisation). 



2. Algebraic and combinatorial preliminaries 
Let k be a field of characteristic zero (in our case it will always be either R or C). 

2.1. Trees. Recall that a tree t is a connected and simply connected graph made out of vertices and edges, the 
sets of which we denote by V(t) and E(t), respectively. 

2.1.1. Planar binary trees. A planar binary tree is a finite oriented tree given an embedding in the plane, such 
that all vertices have exactly two incoming edges and one outgoing edge. An edge can be internal (connecting 
two vertices) or external (with one loose end). The external incoming edges are the leaves. The root edge is 
the unique edge not ending in a vertex. For any planar binary tree t, a partial order on the set of its vertices 
V(t) is defined as follows: u, v e V(t), u < v if and only if there is a path from the root of t through u up to v. 



Y 



The single edge | is the unique planar binary tree without internal vertices. We denote by T 1 "" (resp. T 1 "") 
the set (resp. the linear span) of planar binary trees. A simple grading for such trees is given in terms of the 
number of internal vertices. Above we listed all planar binary trees up to degree three. The number of trees of 

order n is given by the Catalan number c n - , + ", | !n , . The first ones are 1, 1, 2, 5, 14, 42, 132, Alternatively, 

one can use the number of leaves. Observe that for any pair of planar binary trees t\ , fa we can build up a new 
planar binary tree via the grafting operation, t?, := t\ V t%, i.e. by considering the unique Y-shaped planar binary 
tree \, and replacing the left branch (resp. the right branch) by t\ (resp. t-£). 



V 



Y Y V| = V |V Y = V Y V Y : 



Any planar binary tree t ± | obviously expresses itself as t\ V ti in a unique way. The tree t\ (resp. t^) is 
the left part (resp. the right part) of t. The grafting operation V makes T " the free magma algebra with one 
generator: the binary operation V shows no relation of any kind, in particular it is neither commutative nor 
associative. Notice that this product is of degree one with respect to the grading in terms of internal vertices, 
i.e. for two trees t\,t% of degrees n\,ri2, respectively, the product t\ V t% is of degree n\ + n% + 1. However, 
with respect to the leave number grading this product is of degree zero. We call the trees r r , tJ recursively 

a c a u (0) i (0) . (n+1) i,, (») (n+1) (n) ., , 

defined by r r :- \ =: tJ and r r := | V r r , tJ .= r) V |: 



right and left combs, respectively. 
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2.1.2. Planar rooted trees. A planar rooted tree is a finite oriented tree given an embedding in the plane, such 
that all vertices, except one, the root, have arbitrarily many incoming edges and one outgoing edge. The root 
vertex has no outgoing edge, and the leaves have no incoming edges. 



i 



i i V t -i Y 



The single vertex • is the unique rooted tree without edges. Note that we put the root at the bottom of the 
tree. The set (resp. the linear span) of planar non-empty rooted trees will be denoted by T p i (resp. T p i). A 
natural grading for such trees is given in terms of the number of edges. Another one is given by the number of 
vertices. Observe that any rooted tree of degree bigger than zero writes in a unique way: 

t = B + (ti ■ • ■ t n ), 

where B + associates to the forest t\ ■ • • t n the planar tree obtained by grafting all the planar trees tj, j = 1, . . . , n, 
on a common root. 

B + (.)=l 5 + (..)=V, *+(!•) =^, B + (.l)=<}, B + (...)=T- 

Sometimes, one finds the notation t = [t\ •••?„] in the literature 0. Note that the order in which the branch 
trees are displayed has to be taken into account. We introduce the left Butcher product of two planar rooted 
trees t - B + (t\ ■ • -t n ) and u = B + {u\ • ■ • u p ): 

(1) t o^ u := B + (tu\ ■ • ■ Up). 

Hence, it is defined by connecting the root of t via a new edge to the root of u such that t becomes the leftmost 
branch tree. Observe that it is neither associative nor commutative. Moreover, it is clear that any rooted 
tree t = B+(t\ ■ ■ ■ t n ) e T p i of degree bigger than zero uniquely decomposes as t - t\ o^ t 2 . The rooted trees 
recursively defined by ^ 0) := • =: c (0) , and ^ ( " +1) := tfd o-> •, c^' !+1) := • o-» c (,l) are called ladder trees and 
corollas, respectively. 

2.1.3. Knuth's correspondence between planar binary and planar rooted trees. A natural question is how to 
relate the two sets of planar trees just presented. Knuth described in BUI a natural way to do this, known as 
rotation correspondence between planar binary and planar rooted trees. We only give a recursive description 
of this bijection denoted as <1> : T bl " — » T p i, by defining 0(|) := • and: 

(2) «D(fi V t 2 ) := 0(?i) o+ Ofe). 

This map is well-defined and bijective, with its inverse recursively given by: 



(3) ® _1 (* + (fi • • • t n )) = <sr\h) v <sr\B + (t 2 ■ ■ ■ t n )) = 




The first few terms write: 

<D(|) = 



e>(V/) 



o(\y) 



\: o(\\/) - 7 3>(W) = ^ o(\y/) - Y- 



Observe the compatibility with the gradings by the number of internal vertices in T b } " 1 and the number of edges 
in Tpi. This simple bijection implies that the left Butcher product £Q) is also purely magmatic. Left and right 
combs t, , t" map via <1> to the ladder trees {S n) and corollas c ( "\ respectively. For reasons to become clear 
in the sequel, we ask the reader to note the equality between the number of those leaves of a planar binary tree 
which point to the left, and the number of leaves of the corresponding planar rooted tree. 
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2.2. Pre-Lie algebras. Recall that a left pre-Lie algebra (A, >) is a ^-vector space A equipped with an opera- 
tion > : A ® A — > A subject to the following relation: 

(a > b) > c — a > (b > c) - (b > a) > c - b > {a > c). 

See e.g. IT39T1 for a survey on pre-Lie algebras. A genuine example is the pre-Lie algebra of vector fields. Let 
M be a differentiable manifold equipped with a flat, torsion-free connection V. The space of vector fields ^(M) 
can be given the structure of a pre-Lie algebra by defining the product / > g = Vfg. In the case of M = W 
with its canonical flat and torsion-free connection we have that for / = £" =1 fidi and g = YI\=\ Sjdj'- 






d. 



Recall from Chapoton and Livernet |[T3l that the basis of the free pre-Lie algebra in one generator, f{»), can be 
expressed in terms of undecorated, non-planar rooted trees. The set (resp. the linear span) of the latter will be 
denoted by T (resp. 7"). In T too, any rooted tree of degree bigger than zero writes t - B+(t\-- -t n ). However, 
the order in which the branch trees are displayed plays no role. The pre-Lie product in P(») becomes very 
explicit in terms of tree grafting, that is, t\ r\ t% is given by summing over all trees resulting from grafting 
successively the tree t \ to each vertex of t 2 : 

veV(t 2 ) 

where r\ v denotes the grating of the root of t\ via a new edge to vertex v of t%. 



i, »/^v(«/^v«)=i+V. •/^v(«rv(«rv«))=i + 3V+^'+i. 



• r\ • = J , (• r\ •) r\ • 

Any tree can of course be written as a polynomial expression in the generator • using r\ and suitable paren- 
thesizing. See IB El for the description of monomial bases for free pre-Lie algebra. 

Let / be any smooth vector field on R". To go from the free pre-Lie algebra on one generator to the pre-Lie 
algebra of vector fields one applies the elementary differential map Tf : P(») — > ^f(R")» which is defined as the 
unique pre-Lie algebra morphism such that !f/[«] = /. Using standard notations, for t = B + (ti •••?„) e T we 
have: 

r f [t](x) = f n \x){Tf[h}{x), . . .,T f [t n ](x)). 
Elaborating on the introduction of this paper, we mention that Butcher's B-series can be seen as series expan- 
sions b(a) = Z/er ((^(t)/cr(t))t in (the appropriate completion of) f{»). The linear function a on f{») maps trees 
to k, and cr(t) is the symmetry factor of the non-planar tree t. Passing to usual B-series amounts to applying 
the elementary differential map Thf, where hf is a smooth vector field on R" multiplied by the so-called step 
size parameter h: 

B f (a-,y) = Yh v ^T f m). 

Indeed, a canonical B-series consists of a formal power series in the step size parameter h containing elemen- 
tary differentials and arbitrary coefficients encoded in a function a on the set of rooted trees T. Further below 
we will see another natural example of pre-Lie algebra in the context of numerical integration methods. 

As a final remark we mention that composition of B-series leads to what is called Butcher's group, which 
clarified important aspects in the theory of Runge-Kutta methods. Later, Chartier, Hairer and Vilmart intro- 
duced a so-called substitution law for B-series in the context of backward error analysis. See ITT4] l35l for 
concise reviews. 

2.3. Rota-Baxter algebras. Recall that a Rota-Baxter algebra is a ^-algebra A endowed with a ^-linear map 
R : A — > A that satisfies the relation 

(4) R{a)R(b) = R(R(a)b + aR{b) + dab), 

where 9 e k. The map R is called a Rota-Baxter operator of weight 9. The map R :- -GiA - R also is a 
weight 9 Rota-Baxter map. Both images R(A) and R(A) are subalgebras in A. One may think of dU) as a 
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generalized integration by parts identity. Indeed, a simple example is given by the classical integration by 
parts rule showing that the ordinary Riemann integral is a weight zero Rota-Baxter map. Other examples can 
be found for instance in |[T8l[T9"1 . Observe that for an associative Rota-Baxter algebra the two compositions 
a*b '.= R(a)b + aR(b) + Qab and a>b := [R(a), b] - Gba define a new associative product and a pre-Lie product, 
respectively. 

2.4. Dendriform algebras. We introduce the notion of dendriform algebra ||3TI over k, which is a ^-vector 
space D endowed with two bilinear operations -< and > subject to the following three axioms: 

(5) {a < b) < c = a < (b < c + b > c) 

(6) {a > b) < c = a > (b < c) 

(7) a > (b > c) = (a < b + a > b) > c. 

In a commutative dendriform algebra (also known under the name Zinbiel algebra), the left and right operations 
are identified, that is, x > y = y < x. Axioms ©-([7]) imply that for a,b e D the composition: 

(8) a*b:=a<b + a>b 

defines an associative product. Hence, a dendriform algebra is an associative algebra together with a bimodule 
structure on itself, such that the associative product splits into the sum of the left- and right-module struc- 
tures. Moreover, dendriform algebras are at the same time pre-Lie algebras. Indeed, one verifies that the two 
products: 

(9) a> b := a > b — b < a, a<b:=a<b-b>a 
are left pre-Lie and right pre-Lie, respectively. That is, we have: 

(a > b) > c — a > (b > c) - (b > a) > c — b > (a > c), 
(a <b) < c - a <(b < c) - (a < c) < b — a < (c < b). 

These two pre-Lie products vanish if the dendriform algebra is commutative. Recall that a left pre-Lie algebra 
is Lie admissible dEl, that is: 

[a, b] := a < b - b < a 

defines a Lie bracket. An analogous statement holds for right pre-Lie algebras. Moreover, the Lie brackets 
following from the associative operation © and the pre-Lie operations © all define the same Lie bracket. For 
any dendriform algebra A we denote by A = A ffi k.l the corresponding dendriform algebra augmented by a 
unit 1, with the following rules: 

a < I :- a -: I > a I < a :- -: a > I, 

implying a*\ - \*a - a. Note that the equality 1*1 = 1 makes sense, but that 1 < 1 and 1 > 1 are not defined. 

Now suppose that the dendriform A is complete with respect to the topology given by a decreasing filtration 
A - A 1 D A 2 D A 3 D • • • compatible with the dendriform structure, in the sense that A p < A q c A p+q and 
A p > A q c A p+q for any p,q > 1 . In the unital algebra we can then define the exponential and logarithm map 
in terms of the associative product ([8]): 

exp*(.x) :=JV n /n! resp. log*(l + x) := -^(-l)V/n. 

Let L a> (b) :- a > b =: /?>/, (a). Note that Lc^L^ - L at b> and R^R^, = R^ta- We recursively define the set 
of dendriform words in A for fixed elements x\, . . . , x n e A, n £ N by: 



wy(xi,...,x n ) 
wy Ui, . . . , x n ) 

wy(xu . . . , X n ) 



= 1 =: wl\xi x n ) 

= x x <(w { "~ l \x 2 ,...,x n )) 
= (w ( "~\xi,...,x n -i)) > X n . 



In case that x\ = • • • = x n = x we simply write w" (x, . . . , x) = x"' and w" (x, . . . , x) = x" . 
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A simple example of dendriform algebra is given in by the Riemann integral on an algebra f of locally 
integrable functions over R. The dendriform left and right products are: 

if < 8)(0 := fit) If g(s)ds\ (f > g)(t) := U f{s)ds\ g(t). 

With I(f)(t) := L f(s)ds, the dendriform axioms simply encode integration by parts 

if<g)<h = fl(g)l(h)=fl(gl(h)) + fl(l(g)h)=f<(g<h) + f<(g>h) 

(f>g)<h = (I(f)g)I(h) = I(f){gl(h)) = f > (g < h) 

f>ig>h) = I(f)I(g)h = I(fI(g))h + I(I(f)g)h = (f<g)>h + (f>g)>h. 

The pre-Lie product (/ > g) is given by the bracket [/(/), g] . The last example generalizes to any associative 
Rota-Baxter algebra (M, R) of weight 9, giving rise to a dendriform algebra ($[, <, >) defined in terms of: 

a < b := aR(b) + 6ab = -aR(b), a > b := R(a)b. 

The dendriform associative and left pre-Lie products are explicitly given for a, b e 9K by: 

(10) a * b = R(a)b + aR(b) + Bab resp. a > b = [R(a), b] - Gba. 



2.4.1. Dendriform algebra structure on planar trees. In BTl it was shown that planar binary trees different 
from | generate the free dendriform algebra in one generator. The associative product for two trees s = s\ V S2 
and t = h V t% in 7" '" := T bl " — { | } is given recursively by: 

S * t = S[ V (S 2 * t) + * t\) V ? 2 - 

The two terms on the right define the dendriform compositions -< and > respectively. A simple computation 
shows the following link between the dendriform structure and the magmatic product: 

(11) sV t = s> Y <t. 

The tree | can be taken as the unit for the corresponding augmented dendriform algebra. For any dendriform 
algebra A there is a unique morphism F u : T '" — > A. Using (fTTl ). it is recursively given by F a ( \f) - a and: 

F a {t) = F a (t l Vt 2 ) 

= F a {t x >\<t 2 ) 

(12) = F a (h) > a < F a {h). 

This unital dendriform algebra structure on T b J" transfers via Knuth's correspondence between planar binary 
and planar rooted trees (J2j to the dendriform associative product defined for s = s\ o^ s 2 , t = t\ o> ? 2 £ Tpi 
by: 

S • t = S[ 0-^> (S 2 * t) + (S • t[) 0-^> t 2 . 

The tree • is the unit for this product. In the free dendriform algebra with one generator a, the dendriform 
words a" and a" translate into the left and right comb trees, respectively. And via Knuth's correspondence 
they are mapped to the ladder trees and corolla trees, respectively. Observe that for rooted ladder trees we find: 

£ (n)+ € {m) = (f („-l) ^ £(0)) ^ (tfm-l) ^ ^(0)) 

- (t (n ~ l) a, 6 m) ) + {6 n) • 6 m - 1) ) a, 6 0) 



= £ ( • • • ({&-V o, t (m - r) ) o, t W) ) o* • • • ) o* ^ (0) 



)■=() 



For instance: 



:*i = v+i ui = ^ + Y+i U; = lr + 1 
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3. The Magnus expansion 
We start by recalling how to solve the linear initial value problem (IVP): 

4-XQ) = A(t)X(t), X(0) = X , 
at 

with A(t) being, for instance, a matrix valued differentiable function. In fact, let us further simplify the problem 

by assuming scalar- valued functions, i.e. we ignore any commutativity issues. Then the solution of the IVP in 

terms of the exponential map is X(t) = exp( J A(x)cIx)Xq. Expanding the exponential and making use of the 

integration by parts rule immediately yields the identity: 



1 + I A(s\)ds\ + I A{s\) I A{s 2 )ds 2 ds\ + 
Jo Jo Jo 



»( f A(s)ds). 



(13) A(si)dsi + j A(s\) | A{s 2 )ds 2 ds\ +--- = exp| 

The left hand side, known as Dyson-Chen series, corresponds to the formal solution of the integral equation: 

X(t) = l+f A(s)X(s)ds, 



which is associated to the above IVP. This solution is still valid in the case where A is a constant n x n 
matrix. However, this changes drastically in the general non-commutative case. Wilhelm Magnus described in 
a seminal 1954 paper 071 a particular differential equation for the matrix- valued function Q(f; A): 

^n(.s;A) n8« 



Q(s;A) 



-(AM) = A(s) + ^ -lad^ s . A) (A(s)), 

n>0 



e ( a dn(s,A)) 
such that the solution of the IVP writes: 

X(?) = exp(Q(?;A))X . 
It is clear that Q(0; A) = 0, hence Q.(t; A) = J Q(s;A)ds. The B n are the Bernoulli numbers: 

Bo = 1, «i = -r, B 2 = -,B A = -— , ... and B 2k+l = for k > 1. 
2d 3D 

As usual, ad£ (W) stands for the rc-fold iterated Lie bracket [U, [U, • • • [U, W]] •••]. Let us write down the first 
few terms of what is called Magnus' series, Q.(s; A) = X«>o Qi(s',A), following from Picard iteration to solve 
the above recursion: 



Q(s; AA) = AA(s) 



A 3 
+ T 



Jo IJo 



A(x)dx,A(s) 



A(x)dx,A(y) 



dy,A{s) 



A 3 
+ ~2 



\A(x)dx, \A(y)dy,A{s) 
Jo [Jo 



where we introduced a dummy parameter A. Observe that the Lie bracket at order A 2 simply serves to eliminate 
the second term on the righthand side: 

-( I A(s)ds) =2(1 I A( Sl )A(s 2 )ds 2 dsi + I I A(j 2 )A(*i)^ 2 ^i), 
such that up to second order in A: 

exp ( f AA(s) - — \ \A{x)dx,A{s) ds + 0(3) ) = A \ A(s) + A 2 [ \ A{s { )A{s 2 )ds 2 dsi + 0(A 3 ). 
\Jo 2 J [Jo J / Jo Jo Jo 

Let us remark that since its publication, Magnus' paper received much attention and triggered important 
progress in both applied mathematics and physics. We refer the reader to J5] for a comprehensive overview 
of its applications including a concise review of its background. Iserles et al. ETl 1281 1291 used planar tree 
to explore Magnus' expansion. Magnus' result has also been explored from a more algebraic-combinatorial 
perspective using operads, pre-Lie algebras, dendriform algebras and noncommutative symmetric functions 

EdamCElElESl. 
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3.1. The pre-Lie Magnus expansion. Motivated by Spitzer's work P31 in probability theory, Baxter |@), 
see also 0, followed a more general approach to the above IVP suggesting to look at the following fixpoint 
equation: 

(14) X=l+AR(aX), 

in A [[/I]], where A is a commutative unital algebra. The linear map R on A is supposed to satisfy the 
Rota-Baxter relation of scalar weight 6 ©. The solution of ([141 is known as Spitzer's identity. In |[T8l 
we generalized this to a noncommutative Rota-Baxter algebra A by showing that equation (fT4l) is solved by 
X(Aa) = exp(R(fl'(Aa)), where: 

(15) a{Aa) = Aa + Y J ^, (a> {Aa) 

»>0 

denotes the pre-Lie Magnus expansion. Here L x> (y) := x > y where > is the pre-Lie product defined in (fTOT t. 
The prime notation shall remind the reader of Magnus' original differential equation. The first few terms are: 

,1 ,/l 1 „ \ 

Ll (Aa) — Aa — A —a> a + A —(a > a) > a H a > (a > a) + • • • . 

2 \4 12 / 

The reader is refereed to ITT81 [T9l for more details. This leads to the following so-called noncommutative 
Spitzer identity: 

(16) 1 + AR(a) + A 2 R(aR(a)) + A 3 R(aR(aR(a))) + • • • = exp (R{Q.'{Aa))). 

Let us briefly note that proceeding analogously to the example of B-series in Paragraph 12.21 we may abstract 
the pre-Lie Magnus expansion into (the appropriate completion of) free pre-Lie algebra in one generator P(»): 

co'(J3) = Y J (W)/o-(t))t. 

teT 

Now, to go from the free pre-Lie algebra on one generator to the pre-Lie algebra (JR, >) with the pre-Lie product 
defined in (TT0T >. one applies for any a e 9K the map T a '■ P(f) — > &, which is defined as the unique pre-Lie 
algebra morphism such that !F U [•] = a. For rooted trees of higher orders, we must recall that we can choose a 
monomial basis of P(»). Once the monomial basis is fixed, the decomposition of a tree t e T of order n in this 
basis gives rise to a unique polynomial p^(») of order n in the generator •. Then T a \t\ = falp^i*)] '■= pl>( a )> 
and the series co'ifi) in !P(») is mapped to Q.'(a). Of course, this picture demands for more details, such as for 
instance the way to chose the monomial basis. Moreover, we already remarked in lTT8l IT91 that this approach 
will lead to a Magnus expansion with fewer terms than the classical series. A more complete account of this 
perspective on the Magnus expansion will be presented in a forthcoming work. 

3.2. Pre-Lie Magnus expansion and the group of formal flows. Using the dendriform product -< given by 
x < y = xR(y) (we put the weight 6 to zero for simplicity), equation (TToT i in the last paragraph reads: 

1 + R(Aa + A 2 a < a + A 3 a < (a < a) + ■ ■ ■) = exp (R(Q.'(Aa))). 

Using the associative product * =< + > given by x * y = R(x)y + xR(y) and extending the linear map R to the 
associated unital dendriform algebra by setting R(l) := 1, yields: 

X = exp* (D.'(Aa)), 

where: 

X = 2_,(Aa) w = 1 + Aa + A 1 a < a + A 3 a < (a < a) + A 4 a < (a < (a < a)) + • • • 
»>o 
is the solution of the equation X = 1 + Aa < X. 

In |[T8l[T9l we showed in the general dendriform algebra setting the following 



Theorem 1. ( 111811191 ) The element Q! = log* (X(Aa)) in A satisfies the recursive formula similar to (115I) : 

Lq>o , , x vh B„ 



Q , = J^_(^) = V ^l(Aa) 
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with B m the m-th Bernoulli number. The first few terms are: 

,1 ,/l 1 

12 (Aa) - Aa — A —a> a + A I —(a > a) > a H — -a > (a > a) \ + 



3.3. Group of formal flows. In |[20l we adapted a result on the Baker-Campbell-Hausdorff formula by 
Agrachev and Gamkrelidze [fl] to dendriform algebras. Recall that a left pre-Lie algebra A is called chrono- 
logical algebra in yj, and that the left pre-Lie identity rewrites as: 

where as above L clt> : A — » A is defined by L a> (b) = a>b, and where the bracket on the left-hand side is defined 
by [a, b] := a > b - b > a. Recall that as a consequence this bracket satisfies the Jacobi identity. We denote by 
Aue the Lie algebra with the aforementioned bracket. Suppose now that A is a left pre-Lie algebra endowed 
with a decreasing filtration, namely A = A\ D A2 3 A3 d ■ ■ ■ , such that the intersection of the Ay's reduces to 
{0}, and such that A p > A q c A p+q . Suppose moreover that A is complete with respect to this nitration. The 
Baker-Campbell-Hausdorff formula: 

C(a,b) = a +b + -[a,b] + — ([a, [a,b]] + [b, [b, a]]) + ••• 

endows A with a structure of pro-unipotent group. This group admits a more transparent presentation as 
follows. First, introduce a fictitious unit 1 such that 1 > a = a > 1 = a for any a £ A, and define W : A — > A by: 

W(a) := e Lac ~l -l = a + -a>a + -a>(a>a)-\ . 

2 6 

The application W is a bijection. The inverse is the pre-Lie Magnus expansion Q.'(a) introduced in Theorem 
[Q Transferring the BCH product by means of the map W, namely: 

(17) a#b:=w(c(n(a),Q(b))), 

we have W(a)#W(b) = W(C(a,b)) = e L °*e Lb »l - 1, hence W(a)#W(b) - W(a) + e L "»W(b). The product # is 
thus given by the simple formula: 

a#b - a + e La[a> b. 

The inverse is given by a # ~ l = W(-Q.(a)) = e~ La{a)> \ - 1. In particular, when the pre-Lie product > is 
associative, this simplifies to a#b = a > b + a + b and a # ~ l = -j-j^ - 1 = Z«>i( _1 )"«n- If (A, >) and (B, >) are 
two such pre-Lie algebras and if/ : A — > B is a filtration-preserving pre-Lie algebra moiphism, it is immediate 
to check that for any a, b e A we have: 

if/(a#b) = i/<a)#iA(fr). 

In other words, the group of formal flows is a functor from the categoiy of complete filtered pre-Lie algebras to 
the category of groups. Concerning linear dendriform equations (fT8l) in a complete filtered unital dendriform 
algebra A we observe the following interesting fact: for any collection (a\, . . . , a n ) in A with nitration degree 
> 1, the product X = X{a\) * • ■ • * X(a n ) where X{a{) is the solution of the equation X(ai) = 1 + a\ < X{ai) is 
the solution of the linear dendriform equation X = 1 + a < X, with: 



«!#• • -#a n , 



which also writes explicitly: 



»-i 



ci\ + e Ln<a i ,> fl , 2 + \- <? Ln( "i> • ■ • e La(a »-i > a n . 



a closed formula for the magnus expansion 1 1 

4. Linear dendriform equations 

Let us emphasize that, although the series ([731 ) has been known for a long time in the pure pre-Lie context 
Ifn . the results in |[T8l[T9l , interpreting this series as a logarithm, were obtained in the more restricted context 
of dendriform algebras. 

The starting point is the linear dendriform equation in J?l[[/i]] below, where 3\ is now any unital dendriform 
algebra: 

(18) X = l + Aa<X 

for a € {A. Its formal solution is: 

X = \\(Aa)^ = 1 + Aa + A 2 a < a + A 3 a < (a < a) + A 4 a < (a < (a < a)) H . 

n>0 

Remark 2. In this work we deliberately have suppressed any Hopf algebra parlance. However, a few words are 
in order. Recall that the associative algebra (7~ fo ; m , *) is the free algebra generated by the elements |VT, T e T b '" 

ll32l Theorem 3.8]. The subalgebra generated by the right combs t" = \ ( "\ n > is free and is a cocom- 
mutative Hopf subalgebra «H of T b \ n . The left combs t ( " ] = Y"\ n > also belong to <H. It follows that <H 
is isomorphic, as a Hopf algebra, to the Hopf algebra of noncommutative symmetric functions 11221 . In Hopf 
algebraic terms the solution X = 'Z n >o(Aa)^' of (fT8l) is characterized as group-like, that is A(X) = X ® X. 

4.1. A closed form for the logarithm. In this section, discarding the pre-Lie product, we give an explicit 
expression of log*(X) in the free dendriform algebra in one generatoo For this we use Knuth's rotation 
correspondence and consider the representation in terms of planar rooted trees rather than planar binary trees. 
The generator a corresponds to the tree J. Then the solution X is given by the sum of rooted ladder trees, 
X = l + Z„ >0 € {n) =: 1 + L. Here we identify 1 = l^ = .. 



log*(X) = log*(l + L) = ^ 



(-1) 



n+l t *n 



n>0 

«>0 k>0 n+-+i k =n 

i j>0,j=l,...,k 



6 h) *■■■* £ {ik) . 



We have omitted the parameter A. Recall that the degree \ft n ^\ of the tree i^ is equal to n, i.e. its number 
of edges. The series above makes sense in the completion of the free dendriform algebra J{ with respect to 
the grading. Recall that a ladder tree with n edges stands for the dendriform word a < , which in the classical 
example (TTBl corresponds to the n - 1 fold iterated integral a(t)l(al(a ■ ■ ■ I(a)))(t). 

Remark 3. Since X is group-like, the logarithm is given by applying the eulerian idempotent, that is, log*(X) = 
log*(id)(X). 

Theorem 4. The element Q! = log*(X) in J{ is given by the formula: 

£j n £j ( n-l \ ' 

«>0 reT pl \X(t)-1J 

where X(t) denotes the number of leaves of the planar rooted tree t, and |t| its degree, i.e. its number of edges. 
The first few terms are: 



n' 



2 2 3 



= J > v 'I'l't, > T _>-J 'Y + -. 



6 3 6 6 



We recover the q = 1 case of a ^-analog formula by F. Chapoton (see II 1 01 Proposition 5.10]). 
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Proof. First observe that for any ordered composition n = i\ H 1-4 of some positive integer n, the expression 

£fa) • • • • • tS lk) is a linear combination of planar rooted trees with n edges, with coefficients equal to or 1. 
Moreover, for any planar rooted tree r with n edges and k leaves, there is a unique ordered composition 
n = i\ + ■ ■ ■ + 4 of the integer n such that: 

(1) The corresponding *-monomial of ladder trees v- ix > • • • • • v- lk > contains the tree r precisely once. 

(2) For any other composition n - j\ + • • • + j r of the integer n, the •-monomial ft 1 ' • • • • • 1y T > contains 
r with a nontrivial coefficient only if the composition of the integer n is finer than the first one. 

Recall that each leaf of a planar rooted tree is connected to the root by a unique shortest path, i.e. including 
a minimal number of edges. The composition (j\, . . . , 4) is defined as follows: first we number the leaves 
consecutively from left to right. The rightmost leaf is linked down to the root by the path fi^ of length 4 
(i.e. the height of the rightmost leaf). For any s £ {1,. . . ,k — 1}, i s is the length of the path ^ (,) from leaf 
number s to the vertex lying on the unique path joining leaf number s + 1 down to the root. As an example the 
following six-leaved tree with 13 edges: 




is associated to the ordered composition 13 = 2 + 1 + 1+3 + 2 + 4 of its degree. 
Next we write for fixed n > 0: 

k>0 ii+-+i k =n THT pt 

ij>0,j=l,...,k | f | = „ 

The coefficients in the last equality follow by projection, s„(t) := (Z T ,S n ), where Z t (t') = 6 T y. Now, we 
introduce a dummy parameter a in the sum: 

k>0 ii+-+i t =ii 

i'pOJ=l,...,/t 

such that J -j^S n {a)da - S n (-l). Looking first at rooted ladder trees, we observe that in S n (a) only one 
type of rooted ladder tree appears, i.e. the one of length n, and: 

S n (a) = ^ -<** -1 J] &° *■■■* t (k) = J] - a \ k Z - If + non - ladderr °° ted trees - 

i ; >0,j=l,...,/t 

Here we have used that the total number of length k ordered compositions of n is \ k ~_\ Projecting onto the 
ladder part, and integrating, we find the number of rooted ladder trees i^ in the sum S n {-\): 

Sn {& } ) = (z f( „,, s n (-i)> = j x ~ a \ k ~ if a = r ~ (i + a),i " ififa = ^ 

Now consider an arbitrary rooted tree r of degree n with X(r) = k leaves. It is associated with an ordered 
composition n = i\ + • • • + 4 as explained above. Using the fact that considering an ordered composition of n 
finer than the first one is nothing but choosing an ordered composition of i s for any s = 1, . . . ,k, we get: 



«•« -<*.«-'» - f-'-'t- •i/ir-D-^'Cr- 1 .) 



da 



f ^(l+a)' 1 *-"**-*^ 

Jo 
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= 


(- 


-If- 1 

Jo 


0*^(1 


-of 


~ k da 


= 


(- 


-l)*- 1 ^, 


72 — fc + 


1) 




= 


(- 


,!'-*- 


-l)!(n- 


-*)! 






(- 


_l)X(r)-l 









l U(T)-lj 



In the last step we used that k equals the number X(r) of leaves of tree t. u 

Using the inverse of Knuth's correspondence we end up with an equivalent formulation in the planar binary 
tree picture. Recall that a leaf of a planar binary tree is a descent if it is not the leftmost one and if it is pointing 
to the left iflOl . The rotation correspondence yields a bijection between the left-pointing leaves of a planar 
binary tree t and the leaves of its image tree O(f) = r. The leftmost leaf of t (i.e. the one which, by definition, 
is not a descent) is mapped to the leftmost leaf of t. As an immediate consequence we have: 

Corollary 5. In the planar binary tree picture, the element Q! — log*(X) in JK is given by the formula: 



™ ff =Z;Z 



(n-l\ U 
n>0 reT*" \d(t)) 
[j|= n 



where d(t) denotes the number of descents of the planar binary tree t, and \t\ its degree, i.e. its number of 
internal vertices. 

Proof. Follows from the rotation correspondence between planar binary trees and planar rooted trees. □ 



Note that the extra sign (-1)"" 1 in Chapoton's formula ifTOTl could be retrieved starting with the solution 
Y = Y{a) of the equation Y = 1 + Y > a instead of X. Both are linked through the antipode provided the sign 
of the generator is changed, namely Y(a) = S(X(-a)). 

Corollary 6. For any complete filtered dendriform algebra tR = tR\ D J^fe ^ ^3 3 • • ■ and for any a € ZR, the 
element Q'(a) = log* (X(a)) in tA, where X(a) is the solution of the linear dendriform equation X(a) - 1 + a < 
X(a), is given by the formula: 

n>0 re r*;" \d(t)) 

where F a : T' '" — > J{ is the unique dendriform algebra morphism defined in f li2l ), such that F a {\f) = a. 

Hence, using the underlying dendriform algebra structure rather than the pre-Lie one, this result gives a 
closed formula for the Magnus expansion in Theorem [T]in Paragraph 13.21 



5. A FORMULA BY MlELNIK-PLEBANSKl, AND StRICHARTZ 

We would like to give a more precise formulation of Theorem 5] when the dendriform algebra J{ is an 
algebra of matrix-valued functions together with the Riemann integral as a Rota-Baxter operator (of weight 
9 = 0). In the next section we follow (6l [32j. We will see that the continuous Baker-Campbell-Hausdorff 
series, also known as Mielnik-Plebariski-Strichartz formula, gives a closed expression for the classical Magnus 
expansion due to deep structural reasons on the symmetric group algebra, i.e. a natural underlying dendriform 
algebra, which matches with the dendriform algebra structure on planar binary trees. 
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5.1. A dendriform structure on permutations. S„ is the group of permutations of n elements, and k[S n ] its 
group algebra. The concatenation of two permutations cr e S „ and t e S m is the permutation cr x r e S n+m 
obtained by letting cr act on the first n elements and letting t act on the m last elements. An associative product 
is given by: 

cr * t = y w o (cr x r). 
weSh„j„ 

Sh„ m c S n+m stands for the set of (n,m) shuffles. We denote the graded connected algebra with this product 
by <H := ® n > <H„, where «%, - k[S„]. 

The crucial observation is that the product * on permutations splits into * =< + >, where -< and > are 
defined by: 

cr < T - y CO o (cr X t), cr > t = y O) o (cr X t), 

where Sh* m , respectively Sh r ^ m stands for the shuffles co such that co(n + m) = n + m, resp. aj(n) = n + m. It is 
shown in 11331 that < and > endow the augmentation ideal W := (J) ;1>1 *H n with the structure of a dendriform 
algebra. 

Remark: The (dendriform) algebra just described represents only half of the picture. Completing the structure 
amounts to what is known as Malvenuto-Reutenauer Hopf algebra Il38l . which is a graded connected Hopf al- 
gebra on fi. We refrain from giving more details and refer the reader to the references ||6][T7J|23. However, 
let us add that the coproduct splits also into two parts, thus endowing the Hopf algebra fi with the much richer 
structure of bidendriform Hopf algebra, see Foissy's work ETTl . 

5.2. Planar binary trees and permutations. The material presented in this paragraph is mostly borrowed 
from ||32l . See also [6] for a very detailed similar description. A bijective correspondence between permuta- 
tions and planar binary trees with levels is well-known in combinatorics [25]. A planar binary tree with levels 
is a planar binary tree t with, say, n internal vertices together with a bijective decreasing map cp from the poset 
of its internal vertices into {1, . . . , n\. Such a tree admits a graphical realization by drawing the internal vertices 
at the prescribed levels, with level 1 being the top one and level n being the deepest one. Any planar binary tree 
with levels (t, ip) gives rise to two such trees (t\,<p\) and (?2, yi), where t = t\ V ?2 and tpi is the "standardized" 
restriction of the injection tp to the internal vertices of ?,-, namely its composition on the left with the unique 
increasing bijection from its image onto {1, . . . , |?;|}, i - 1,2. 

To any such tree (t, tp) we can associate a permutation cr t>ip as follows: cr t ^(i) is the level of the internal 
vertex Uj situated between leaves /,• and Z r+ i (the leftmost being the first and the rightmost being number n + 1). 
This correspondence P is a bijection, the inverse of which is recursively given as follows: the permutation 
cr e S n gives rise to two sequences of integers: the sequence before n and the sequence after n in (<ri, . . . , cr„). 
One of them may be empty. By "standardizing" the integers in each sequence, they form a permutation. For 
instance (341625) gives the two sequences (341) and (25), which, after standardizing, give (231) and (12). By 
induction these two permutations give rise to two trees with levels. The grafting V of the two trees (in the order 
given above) gives the underlying tree of P~ l (cr), and the original permutation is used to determine the levels 
of each vertex, namely (p(Uj) = cr~ l (j). 

Recall that the descent set of a permutation cr e S„ is the subset D(cr) c {1, . . . , n - 1 } of indices i such that 
cr(i) > cr(i +1). We denote by d(cr) the cardinality of the descent set. The following lemma, the proof of which 
is simple, establishes a link with descents of planar binary trees: 

Lemma 7. For any planar binary tree with levels (t, ip), the correspondence 

j h-> lj + \ : { 1, . . . , n\ — > {leaves of t) 

restricts to a bijection from the descent set D(cr ttf ) onto the descent set oft. 

Forgetting the levels, the bijection P~ l provides a surjective map ifr„ : S n — > [T b ™) for any n > 0. As a 
corollary of Lemma [7] we have: 

(22) d(cr) = d(tfr(cr)) 
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for any cr e S n . Dually, we get a linear injection i]/* n : [T b J n ) — > k[S n ] given by: 

iA„(o-)=f 

These maps together give a degree zero linear injection if/* from T b '" to the dendriform algebra ( H. The 
following theorem is due to Loday and Ronco ( 11321 Theorem 3.1], IT331 Proposition 5.3]): 



Theorem 8. [33] The linear injection if/* defined above is a unital dendriform algebra morphism. 
Remark 9. The map if/* is even a bidendriform Hopf algebra moiphism lETI . 

5.3. The Mielnik-Plebaiiski-Strichartz formula. Let Si be an algebra of locally integrable functions of 
one real variable, with values in some topological not necessarily commutative algebra, for example Mn(C), 
N G N. The weight zero Rota-Baxter map R : 3\ — > Si given by Ra{s) := L a{u)du endows Si with a 
dendriform algebra structure, as explained in Paragraph 12. 41 above. Fix some a in Si, suppose s > 0, and define 
Fcr{a) € Si for any permutation cr e S ,„ n > 1 as follows: 

F a (o~)(s) := — /■■■I aiua-i) ■ ■ ■ a(uo- n ) du\ ■ ■ ■ du n . 

0<«„<--<Ul<.V 

In particular, F fl ((l)) = a, F a ((12)) - a/?(a), F a ((21)) = R{a)a, F fl ((123)) = aR(aR(a)) and F a ((321)) = 
R(R(a)a)a. Obsei"ve that there are permutations cr e S„, n > 2, which, strictly speaking, can not be written as 
iterated operators, e.g. (231) and (132). However, they can be written as iterated integrals in the above sense. 
Therefore, the following only applies to the special case of the dendriform algebra structure defined in terms 
of Ra(s) := J a(u)du on Si. The crucial point is to transfer the dendriform calculations to the simplicies 

forming the integration domains. The correspondence F a obviously extends uniquely to a linear map from the 
augmentation ideal of <H to Si. 

Before stating the theorem, let us introduce some notations: for any positive integer ,/V we denote by AjL c 

[0, s\ N the simplex {u = (ui Un), < Un < • • • < u\ < s}. The symmetric group S # acts on [0, s] N by 

permutation of the variables, namely: 

cr.(ui,...,u N ) := (u a .-i,...,u (r -i). 

We have then: 

CT.A S N = {U = (ll\, . . . , Un), < Ua- N < ■ ■ ■ < U (Tl < s}. 

Theorem 10. F a is a dendriform algebra morphism from the augmentation ideal of 1-1 to Si, linked to the 
dendriform algebra morphism F a of Corollary\6\by: 

F a = F a oif/*. 
Proof. By direct computation: take cr e S„ and t e S ,„. Then, 

R(F a (cr<T))(s) = J] R(F a (coo(crXT))(s) 

a>eSh~(n,m) 



2j I '] a ( u <r\ ) " ' a ( u o-„ )"(m»+ti ) • • • a(u n+Tm ) du 

in (n,m) %+ m w ! 

2] I • ■ ■ J afro-! ) • • • a{u (Tn )a(u n+ n ) • • • a(««+r„, ) <^« 



weSh (n,ni) 



\-\ a(u ai )--- a(Uo- n )a(u n+Tl )■■■ a(u n+Tm ) du 

J J [J &i _1 .A„ + „, 



weSh 2 (n,m) 
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a(u (Tl ) • • • a(u (Tn )a(v Tl ) • • • a{v Tm )dudv 



J ' ' J°< u » 



<--<m<s 



= R(F a (cr)R(F a (T))) 

= R(F a (cr) < F a {r)), 

hence F a (cr < t) = F a {cr) < F a (r). The corresponding computation for > is entirely similar and left to the 
reader, using Sh'(«, m) instead of Sh 2 (n, m). The second statement is an immediate consequence of Theorem 
[8]and the freeness of the dendriform algebra of planar binary trees. □ 

Corollary 11. The element Q.'(a) - log* (X(a)) in tft, where X(a) is the solution of the linear dendriform 
equation X(a) = 1 + a < X(a), is formally given by the series: 



(23) «'=>.- > , ^T^T^Foio-) 



1 (^ 

" L^nLi tn-X\ **• 



(-l)d(o-) 

_i (n-\\ 
h>0 <T£S„ \d(cr)) 



Proof This is an immediate consequence of Corollary l6l Theorem [101 and Eq (l22l). □ 

Corollary 12 (Mielnik-Plebariski-Strichartz formula 02] @6]). The Magnus element Cl(a) = R(Cl'(a)) is 
formally given by the series: 

(24) &(a)(s) = 2_j2l (n-\\ I " ' I a( - u ^ ^ " ' a ( u o-J du i ■ ■ ■ du n- 

Remark : As Q(o r ) is a Lie element, we can use the Dynkin-Specht-Wever theorem, so that we recover the 
formula in its original setting: 

(_-\\d(cr) p r 

Q(a)(s) - ^ ^ 2 (n-i\ \"\ W^i)'t fl (^)""[ fl K-p a("(r„)]---]]^«i •■•rf« n . 

The order of the m/s is reversed compared to the original, since Z = exp(Q.(a)) solves the initial value problem 
Z = aZ rather than Z = Za as in Roll . 
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